demo daphnia instructions
- Load the R scripts in chunk #1
- Run the chunk with the function demo_daphnia(). Run it from the console and not from the notebook otherwise it will take longer will not be able to see the animation.
# chunk #1
library(ReacTran)
library(simecol)
source("../R/daphnia_deb_phyto.R")
source("../R/daphnia_deb_equations.R")
source("../R/daphnia_deb_lifeequations.R")
source("../R/daphnia_deb.R")
source("../R/daphnia_deb_ibm.R")
source("../R/deb_ibm.R")
time= 1 ref. Vol= 10 abu= 10
time= 2 ref. Vol= 20 abu= 9.2
time= 3 ref. Vol= 20 abu= 8.85
time= 4 ref. Vol= 20 abu= 8.1
time= 5 ref. Vol= 20 abu= 7.8
time= 6 ref. Vol= 2.083333 abu= 49.92
time= 7 ref. Vol= 4.166667 abu= 47.04
time= 8 ref. Vol= 4.166667 abu= 44.64
time= 9 ref. Vol= 4.166667 abu= 42
time= 10 ref. Vol= 0.6127451 abu= 168.096
time= 11 ref. Vol= 1.22549 abu= 153.408
time= 12 ref. Vol= 1.22549 abu= 144.432
time= 13 ref. Vol= 1.22549 abu= 136.272
time= 14 ref. Vol= 0.4589851 abu= 235.3018
time= 15 ref. Vol= 0.9179702 abu= 215.6933
time= 16 ref. Vol= 0.9179702 abu= 162.3146
time= 17 ref. Vol= 0.9179702 abu= 124.187
time= 18 ref. Vol= 0.9179702 abu= 108.936
time= 19 ref. Vol= 2.753911 abu= 70.8084
time= 20 ref. Vol= 2.753911 abu= 39.21696
time= 21 ref. Vol= 5.507821 abu= 26.14464
time= 22 ref. Vol= 11.01564 abu= 14.16168
time= 23 ref. Vol= 22.03128 abu= 8.98722
time= 24 ref. Vol= 22.03128 abu= 5.08368
time= 25 ref. Vol= 66.09385 abu= 2.85957
time= 26 ref. Vol= 66.09385 abu= 1.57352
time= 27 ref. Vol= 198.2816 abu= 0.93806
time= 28 ref. Vol= 198.2816 abu= 0.57494
time= 29 ref. Vol= 396.5631 abu= 0.39338
time= 30 ref. Vol= 793.1262 abu= 0.2219067
time= 31 ref. Vol= 793.1262 abu= 0.1386917
time= 32 ref. Vol= 1586.252 abu= 0.08573667
time= 33 ref. Vol= 3172.505 abu= 0.05043333
time= 34 ref. Vol= 3172.505 abu= 0.04854208
time= 35 ref. Vol= 3172.505 abu= 0.04665083
time= 36 ref. Vol= 3172.505 abu= 0.04539
time= 37 ref. Vol= 3172.505 abu= 0.04318354
time= 38 ref. Vol= 3172.505 abu= 0.04129229
time= 39 ref. Vol= 3172.505 abu= 0.03877062
time= 40 ref. Vol= 3172.505 abu= 0.03687937
time= 41 ref. Vol= 3172.505 abu= 0.03435771
time= 42 ref. Vol= 3172.505 abu= 0.03246646
time= 43 ref. Vol= 6345.01 abu= 0.03026
time= 44 ref. Vol= 6345.01 abu= 0.02821115
time= 45 ref. Vol= 6345.01 abu= 0.02679271
time= 46 ref. Vol= 6345.01 abu= 0.02553187
time= 47 ref. Vol= 369.7558 abu= 0.2731532
time= 48 ref. Vol= 739.5116 abu= 0.2569263
time= 49 ref. Vol= 739.5116 abu= 0.2474606
time= 50 ref. Vol= 739.5116 abu= 0.2366427
time= 51 ref. Vol= 331.6196 abu= 0.2834573
time= 52 ref. Vol= 663.2391 abu= 0.2593333
time= 53 ref. Vol= 663.2391 abu= 0.2457635
time= 54 ref. Vol= 663.2391 abu= 0.2352093
time= 55 ref. Vol= 292.1758 abu= 0.3627952
time= 56 ref. Vol= 292.1758 abu= 0.3525274
time= 57 ref. Vol= 584.3517 abu= 0.3080337
time= 58 ref. Vol= 584.3517 abu= 0.3028998
time= 59 ref. Vol= 584.3517 abu= 0.2874981
time= 60 ref. Vol= 62.23127 abu= 1.783669
time= 61 ref. Vol= 62.23127 abu= 1.671185
time= 62 ref. Vol= 62.23127 abu= 1.606909
time= 63 ref. Vol= 124.4625 abu= 1.494425
time= 64 ref. Vol= 27.47518 abu= 3.748839
time= 65 ref. Vol= 54.95035 abu= 3.566856
time= 66 ref. Vol= 54.95035 abu= 3.403072
time= 67 ref. Vol= 54.95035 abu= 3.239288
time= 68 ref. Vol= 19.55529 abu= 5.420529
time= 69 ref. Vol= 19.55529 abu= 5.164844
time= 70 ref. Vol= 39.11057 abu= 4.909159
time= 71 ref. Vol= 39.11057 abu= 4.602336
time= 72 ref. Vol= 15.21812 abu= 6.636825
time= 73 ref. Vol= 6.503471 abu= 14.91511
time= 74 ref. Vol= 13.00694 abu= 14.60759
time= 75 ref. Vol= 13.00694 abu= 14.14629
time= 76 ref. Vol= 5.993982 abu= 16.1829
time= 77 ref. Vol= 2.814076 abu= 31.98208
time= 78 ref. Vol= 5.628152 abu= 29.84994
time= 79 ref. Vol= 5.628152 abu= 28.78387
time= 80 ref. Vol= 5.628152 abu= 27.36244
time= 81 ref. Vol= 1.803895 abu= 65.41401
time= 82 ref. Vol= 1.803895 abu= 63.19659
time= 83 ref. Vol= 1.803895 abu= 59.31609
time= 84 ref. Vol= 1.803895 abu= 58.20738
time= 85 ref. Vol= 1.803895 abu= 95.34923
time= 86 ref. Vol= 0.5352804 abu= 177.4771
time= 87 ref. Vol= 1.070561 abu= 162.5316
time= 88 ref. Vol= 1.070561 abu= 154.1248
time= 89 ref. Vol= 1.070561 abu= 144.7839
time= 90 ref. Vol= 1.070561 abu= 132.6408
time= 91 ref. Vol= 2.141122 abu= 88.73854
time= 92 ref. Vol= 2.141122 abu= 64.91925
time= 93 ref. Vol= 4.282243 abu= 41.567
time= 94 ref. Vol= 4.282243 abu= 29.89088
time= 95 ref. Vol= 8.564486 abu= 19.14884
time= 96 ref. Vol= 8.564486 abu= 12.25993
time= 97 ref. Vol= 25.69346 abu= 7.472719
time= 98 ref. Vol= 25.69346 abu= 4.787211
time= 99 ref. Vol= 51.38692 abu= 2.88011
time= 100 ref. Vol= 102.7738 abu= 1.517896
time= 101 ref. Vol= 205.5477 abu= 0.9243598
time= 102 ref. Vol= 205.5477 abu= 0.5302906
time= 103 ref. Vol= 616.643 abu= 0.2773079
time= 104 ref. Vol= 616.643 abu= 0.1670334
time= 105 ref. Vol= 1233.286 abu= 0.15406
time= 106 ref. Vol= 1233.286 abu= 0.143519
time= 107 ref. Vol= 1233.286 abu= 0.1362214
time= 108 ref. Vol= 1233.286 abu= 0.1329781
time= 109 ref. Vol= 1233.286 abu= 0.1248697
time= 110 ref. Vol= 1233.286 abu= 0.1183829
time= 111 ref. Vol= 1233.286 abu= 0.1151396
time= 112 ref. Vol= 1233.286 abu= 0.1086528
time= 113 ref. Vol= 1233.286 abu= 0.09811187
time= 114 ref. Vol= 1233.286 abu= 0.09081429
time= 115 ref. Vol= 1233.286 abu= 0.08594924
time= 116 ref. Vol= 2466.572 abu= 0.08027335
time= 117 ref. Vol= 1209.104 abu= 0.09097646
time= 118 ref. Vol= 145.325 abu= 0.6261827
time= 119 ref. Vol= 290.65 abu= 0.6124205
time= 120 ref. Vol= 290.65 abu= 0.5780148
time= 121 ref. Vol= 290.65 abu= 0.5573714
time= 122 ref. Vol= 74.90979 abu= 1.428385
time= 123 ref. Vol= 74.90979 abu= 1.334939
time= 124 ref. Vol= 149.8196 abu= 1.148048
time= 125 ref. Vol= 149.8196 abu= 1.101325
time= 126 ref. Vol= 43.0516 abu= 2.438934
time= 127 ref. Vol= 43.0516 abu= 2.346022
time= 128 ref. Vol= 86.10321 abu= 2.229882
time= 129 ref. Vol= 86.10321 abu= 2.171812
time= 130 ref. Vol= 30.00112 abu= 3.233213
time= 131 ref. Vol= 10.38101 abu= 9.632974
time= 132 ref. Vol= 20.76202 abu= 9.054996
time= 133 ref. Vol= 20.76202 abu= 8.814172
time= 134 ref. Vol= 20.76202 abu= 8.332523
time= 135 ref. Vol= 5.832028 abu= 17.6611
time= 136 ref. Vol= 11.66406 abu= 16.97523
time= 137 ref. Vol= 11.66406 abu= 16.46083
time= 138 ref. Vol= 11.66406 abu= 15.60349
time= 139 ref. Vol= 4.21085 abu= 22.56076
time= 140 ref. Vol= 8.421701 abu= 21.37335
time= 141 ref. Vol= 8.421701 abu= 20.66091
time= 142 ref. Vol= 8.421701 abu= 19.71098
time= 143 ref. Vol= 2.089752 abu= 47.37404
time= 144 ref. Vol= 1.029435 abu= 90.34086
time= 145 ref. Vol= 2.058869 abu= 87.42664
time= 146 ref. Vol= 2.058869 abu= 82.5696
time= 147 ref. Vol= 0.8990695 abu= 119.0119
time= 148 ref. Vol= 0.8990695 abu= 189.0844
time= 149 ref. Vol= 0.8990695 abu= 180.1863
time= 150 ref. Vol= 0.8990695 abu= 173.5127
time= 151 ref. Vol= 0.8990695 abu= 161.2779
time= 152 ref. Vol= 0.8990695 abu= 170.176
time= 153 ref. Vol= 0.8990695 abu= 114.5629
time= 154 ref. Vol= 1.798139 abu= 87.86863
time= 155 ref. Vol= 1.798139 abu= 68.40406
time= 156 ref. Vol= 1.798139 abu= 58.39371
time= 157 ref. Vol= 5.394417 abu= 33.36783
time= 158 ref. Vol= 5.394417 abu= 19.64995
time= 159 ref. Vol= 16.18325 abu= 11.86412
time= 160 ref. Vol= 16.18325 abu= 7.476866
time= 161 ref. Vol= 48.54975 abu= 4.078291
time= 162 ref. Vol= 48.54975 abu= 2.451094
time= 163 ref. Vol= 97.0995 abu= 1.503612
time= 164 ref. Vol= 194.199 abu= 0.8341958
time= 165 ref. Vol= 194.199 abu= 0.5200851
time= 166 ref. Vol= 582.597 abu= 0.288364
time= 167 ref. Vol= 1165.194 abu= 0.1287339
time= 168 ref. Vol= 2330.388 abu= 0.07981503
time= 169 ref. Vol= 2330.388 abu= 0.06865809
time= 170 ref. Vol= 2330.388 abu= 0.0622214
time= 171 ref. Vol= 2330.388 abu= 0.05964672
time= 172 ref. Vol= 2330.388 abu= 0.05664293
time= 173 ref. Vol= 2330.388 abu= 0.05492648
time= 174 ref. Vol= 2330.388 abu= 0.05321002
time= 175 ref. Vol= 2330.388 abu= 0.04977712
time= 176 ref. Vol= 2330.388 abu= 0.0523518
time= 177 ref. Vol= 2330.388 abu= 0.04891889
time= 178 ref. Vol= 2330.388 abu= 0.04763155
time= 179 ref. Vol= 2330.388 abu= 0.04720244
time= 180 ref. Vol= 2330.388 abu= 0.06822898
time= 181 ref. Vol= 1049.724 abu= 0.09335784
time= 182 ref. Vol= 294.0404 abu= 0.350292
time= 183 ref. Vol= 588.0809 abu= 0.3298866
time= 184 ref. Vol= 268.5301 abu= 0.3649498
time= 185 ref. Vol= 537.0601 abu= 0.3426059
time= 186 ref. Vol= 144.7601 abu= 0.6009941
time= 187 ref. Vol= 289.5203 abu= 0.5802702
time= 188 ref. Vol= 289.5203 abu= 0.5630003
time= 189 ref. Vol= 289.5203 abu= 0.5422763
time= 190 ref. Vol= 121.6472 abu= 0.8878134
time= 191 ref. Vol= 121.6472 abu= 0.8384904
time= 192 ref. Vol= 243.2944 abu= 0.780947
time= 193 ref. Vol= 243.2944 abu= 0.7357343
time= 194 ref. Vol= 86.89084 abu= 1.242939
time= 195 ref. Vol= 36.6628 abu= 2.863938
time= 196 ref. Vol= 73.32561 abu= 2.673009
time= 197 ref. Vol= 73.32561 abu= 2.522993
time= 198 ref. Vol= 26.28158 abu= 3.5386
time= 199 ref. Vol= 26.28158 abu= 5.441073
time= 200 ref. Vol= 26.28158 abu= 5.288875
# chunk #2
demo_daphnia <- function() {
## load individual-level model (DEB and phyto, differential equations)
submodel <- daphnia_deb_phyto()
## test one outer step alone
m <- sim(submodel)
out(m)
### load individual-based model (IBM, discrete, population dynamics)
deb_ibm <- daphnia_deb_ibm()
initfunc(deb_ibm) <- function(obj) {
p <- as.list(parms(obj))
parms(obj)["WAM"] <- length2weight(p$SAM, p$l2w) # weight at maturity
parms(obj)["WON"] <- length2weight(p$SON, p$l2w) # weight of neonates
age <- sample(1:30, 100, repl=TRUE)
len <- p$SON + 2 * age/30 ## simple linear approximation
weight <- length2weight(len, p$l2w)
initlist <- list(X = c(x1=0, x2=0, x3=0.25),
Z = data.frame(
age = age,
weight = weight,
egg = sample(1:5, 100, repl=TRUE),
eggage = runif(100, min=0, max=5),
length = len,
eggstorage = 0
),
refVol = 10 # 10L
)
init(obj) <- initlist #list2vec(initlist)
obj
}
observer(deb_ibm) <- function (statelist) {
X <- statelist$X
Z <- statelist$Z
refVol <- statelist$refVol
hist(Z$length, col="red", breaks=seq(0.5,3,0.25), ylim=c(0,150))
c(X["x1"], X["x2"], X["x3"], abundance = nrow(Z)/refVol,
age = mean(Z$age), weight = mean(Z$weight), eggs = mean(Z$egg),
length = mean(Z$length), refVol = refVol)
}
#deb_ibm <- initialize(deb_ibm)
times(deb_ibm)["to"] <- 200
parms(deb_ibm)$samplesize <- c(100, 200)
parms(deb_ibm)$life.span <- 60
deb_ibm <- sim(deb_ibm)
o <- out(deb_ibm)
### plot results
par(mfrow=c(3,1))
plot(o$time, o$x1, type="l",
main="Phytoplankton", col="red", ylim=c(0,0.5),xlab="Day",ylab="mg C / L")
lines(o$time, o$x2, col="blue")
lines(o$time, o$x3, col="green")
plot(o$time, o$abundance, type="l", main="Abundance",xlab="Day", ylab="Ind / L")
bm.mgc <- o$weight * o$abund / 1000
plot(o$time, bm.mgc , type="l", main="Biomass",xlab="Day", ylab="mg C / L")
}
demo_daphnia()
time= 1 ref. Vol= 10 abu= 10
time= 2 ref. Vol= 20 abu= 9.8
time= 3 ref. Vol= 20 abu= 9.3
time= 4 ref. Vol= 20 abu= 8.8
time= 5 ref. Vol= 20 abu= 8.1
time= 6 ref. Vol= 2.053388 abu= 49.187
time= 7 ref. Vol= 4.106776 abu= 46.752
time= 8 ref. Vol= 1.740159 abu= 59.76464
time= 9 ref. Vol= 3.480319 abu= 56.31668
time= 10 ref. Vol= 0.6481041 abu= 160.4681
time= 11 ref. Vol= 0.6481041 abu= 154.2962
time= 12 ref. Vol= 0.6481041 abu= 171.2688
time= 13 ref. Vol= 0.6481041 abu= 160.4681
time= 14 ref. Vol= 0.6481041 abu= 197.4991
time= 15 ref. Vol= 0.6481041 abu= 185.1555
time= 16 ref. Vol= 1.296208 abu= 138.8666
time= 17 ref. Vol= 1.296208 abu= 101.8355
time= 18 ref. Vol= 2.592416 abu= 74.83366
time= 19 ref. Vol= 2.592416 abu= 46.6746
time= 20 ref. Vol= 5.184832 abu= 30.4735
time= 21 ref. Vol= 5.184832 abu= 21.02286
time= 22 ref. Vol= 15.5545 abu= 12.53657
time= 23 ref. Vol= 15.5545 abu= 7.39336
time= 24 ref. Vol= 31.10899 abu= 4.436016
time= 25 ref. Vol= 62.21799 abu= 2.507313
time= 26 ref. Vol= 124.436 abu= 1.52689
time= 27 ref. Vol= 124.436 abu= 0.892025
time= 28 ref. Vol= 373.3079 abu= 0.5303932
time= 29 ref. Vol= 373.3079 abu= 0.3402017
time= 30 ref. Vol= 746.6159 abu= 0.206264
time= 31 ref. Vol= 746.6159 abu= 0.1366164
time= 32 ref. Vol= 1493.232 abu= 0.1178652
time= 33 ref. Vol= 1493.232 abu= 0.1111683
time= 34 ref. Vol= 1493.232 abu= 0.1051411
time= 35 ref. Vol= 1493.232 abu= 0.0984442
time= 36 ref. Vol= 1493.232 abu= 0.09040794
time= 37 ref. Vol= 1493.232 abu= 0.08438074
time= 38 ref. Vol= 1493.232 abu= 0.08371105
time= 39 ref. Vol= 1493.232 abu= 0.07567479
time= 40 ref. Vol= 1493.232 abu= 0.06830822
time= 41 ref. Vol= 2986.464 abu= 0.0636204
time= 42 ref. Vol= 2986.464 abu= 0.06161133
time= 43 ref. Vol= 669.6107 abu= 0.1403801
time= 44 ref. Vol= 1339.221 abu= 0.1344065
time= 45 ref. Vol= 195.5068 abu= 0.5114913
time= 46 ref. Vol= 391.0135 abu= 0.4859167
time= 47 ref. Vol= 186.1969 abu= 0.5639191
time= 48 ref. Vol= 186.1969 abu= 0.5370658
time= 49 ref. Vol= 80.25729 abu= 1.183693
time= 50 ref. Vol= 160.5146 abu= 1.146313
time= 51 ref. Vol= 77.9197 abu= 1.360375
time= 52 ref. Vol= 77.9197 abu= 1.334707
time= 53 ref. Vol= 77.9197 abu= 1.732553
time= 54 ref. Vol= 77.9197 abu= 1.668384
time= 55 ref. Vol= 77.9197 abu= 1.540047
time= 56 ref. Vol= 77.9197 abu= 2.027729
time= 57 ref. Vol= 77.9197 abu= 2.41274
time= 58 ref. Vol= 27.92821 abu= 3.75964
time= 59 ref. Vol= 27.92821 abu= 3.652221
time= 60 ref. Vol= 27.92821 abu= 4.798016
time= 61 ref. Vol= 27.92821 abu= 4.511568
time= 62 ref. Vol= 10.01011 abu= 9.9899
time= 63 ref. Vol= 20.02022 abu= 9.790102
time= 64 ref. Vol= 8.483144 abu= 12.14172
time= 65 ref. Vol= 8.483144 abu= 11.78808
time= 66 ref. Vol= 8.483144 abu= 21.80795
time= 67 ref. Vol= 8.483144 abu= 20.62914
time= 68 ref. Vol= 8.483144 abu= 23.57616
time= 69 ref. Vol= 3.70443 abu= 28.0745
time= 70 ref. Vol= 3.70443 abu= 39.41227
time= 71 ref. Vol= 3.70443 abu= 41.3019
time= 72 ref. Vol= 3.70443 abu= 40.22211
time= 73 ref. Vol= 3.70443 abu= 48.86042
time= 74 ref. Vol= 1.722991 abu= 56.87785
time= 75 ref. Vol= 1.722991 abu= 71.96789
time= 76 ref. Vol= 1.722991 abu= 67.90519
time= 77 ref. Vol= 1.722991 abu= 88.2187
time= 78 ref. Vol= 1.722991 abu= 84.73639
time= 79 ref. Vol= 0.640517 abu= 131.1441
time= 80 ref. Vol= 1.281034 abu= 129.5828
time= 81 ref. Vol= 1.281034 abu= 150.6596
time= 82 ref. Vol= 1.281034 abu= 124.8991
time= 83 ref. Vol= 1.281034 abu= 107.7255
time= 84 ref. Vol= 2.562068 abu= 67.13327
time= 85 ref. Vol= 2.562068 abu= 46.44686
time= 86 ref. Vol= 5.124136 abu= 30.44416
time= 87 ref. Vol= 5.124136 abu= 21.07673
time= 88 ref. Vol= 15.37241 abu= 10.14805
time= 89 ref. Vol= 30.74481 abu= 5.854646
time= 90 ref. Vol= 30.74481 abu= 3.480262
time= 91 ref. Vol= 61.48963 abu= 2.276807
time= 92 ref. Vol= 122.9793 abu= 1.39861
time= 93 ref. Vol= 122.9793 abu= 0.8375396
time= 94 ref. Vol= 368.9378 abu= 0.4960186
time= 95 ref. Vol= 368.9378 abu= 0.3903097
time= 96 ref. Vol= 368.9378 abu= 0.3659154
time= 97 ref. Vol= 368.9378 abu= 0.3577839
time= 98 ref. Vol= 368.9378 abu= 0.3442315
time= 99 ref. Vol= 368.9378 abu= 0.3333896
time= 100 ref. Vol= 368.9378 abu= 0.3062847
time= 101 ref. Vol= 368.9378 abu= 0.2981533
time= 102 ref. Vol= 368.9378 abu= 0.2981533
time= 103 ref. Vol= 368.9378 abu= 0.2873113
time= 104 ref. Vol= 368.9378 abu= 0.2710484
time= 105 ref. Vol= 368.9378 abu= 0.3875992
time= 106 ref. Vol= 368.9378 abu= 0.4445194
time= 107 ref. Vol= 368.9378 abu= 0.4499404
time= 108 ref. Vol= 368.9378 abu= 0.425546
time= 109 ref. Vol= 41.54705 abu= 2.238426
time= 110 ref. Vol= 83.09409 abu= 2.118081
time= 111 ref. Vol= 40.5337 abu= 2.368399
time= 112 ref. Vol= 81.06741 abu= 2.09702
time= 113 ref. Vol= 25.57331 abu= 4.692391
time= 114 ref. Vol= 25.57331 abu= 4.457772
time= 115 ref. Vol= 25.57331 abu= 5.083424
time= 116 ref. Vol= 25.57331 abu= 4.770598
time= 117 ref. Vol= 25.57331 abu= 7.625136
time= 118 ref. Vol= 25.57331 abu= 7.507826
time= 119 ref. Vol= 11.16739 abu= 9.491926
time= 120 ref. Vol= 11.16739 abu= 9.13374
time= 121 ref. Vol= 22.33477 abu= 8.417368
time= 122 ref. Vol= 4.370797 abu= 24.70945
time= 123 ref. Vol= 4.370797 abu= 24.70945
time= 124 ref. Vol= 4.370797 abu= 25.85341
time= 125 ref. Vol= 4.370797 abu= 24.93824
time= 126 ref. Vol= 1.783999 abu= 56.61439
time= 127 ref. Vol= 3.567998 abu= 53.8117
time= 128 ref. Vol= 1.585777 abu= 61.79937
time= 129 ref. Vol= 3.171553 abu= 61.16876
time= 130 ref. Vol= 1.026393 abu= 110.0943
time= 131 ref. Vol= 1.026393 abu= 100.3515
time= 132 ref. Vol= 2.052785 abu= 96.45431
time= 133 ref. Vol= 2.052785 abu= 91.09574
time= 134 ref. Vol= 0.986916 abu= 95.2462
time= 135 ref. Vol= 0.986916 abu= 132.7367
time= 136 ref. Vol= 0.986916 abu= 112.4716
time= 137 ref. Vol= 1.973832 abu= 85.11362
time= 138 ref. Vol= 1.973832 abu= 70.42139
time= 139 ref. Vol= 1.973832 abu= 54.20927
time= 140 ref. Vol= 5.921496 abu= 27.86458
time= 141 ref. Vol= 5.921496 abu= 17.0565
time= 142 ref. Vol= 11.84299 abu= 11.31471
time= 143 ref. Vol= 23.68598 abu= 6.332859
time= 144 ref. Vol= 47.37197 abu= 4.137468
time= 145 ref. Vol= 47.37197 abu= 2.70202
time= 146 ref. Vol= 94.74394 abu= 1.519886
time= 147 ref. Vol= 189.4879 abu= 0.9499288
time= 148 ref. Vol= 189.4879 abu= 0.8813229
time= 149 ref. Vol= 189.4879 abu= 0.8602133
time= 150 ref. Vol= 189.4879 abu= 0.8391038
time= 151 ref. Vol= 189.4879 abu= 0.7599431
time= 152 ref. Vol= 189.4879 abu= 0.717724
time= 153 ref. Vol= 189.4879 abu= 0.6649502
time= 154 ref. Vol= 189.4879 abu= 0.6385633
time= 155 ref. Vol= 189.4879 abu= 0.6121764
time= 156 ref. Vol= 189.4879 abu= 0.5857895
time= 157 ref. Vol= 189.4879 abu= 0.5277382
time= 158 ref. Vol= 189.4879 abu= 0.8496586
time= 159 ref. Vol= 189.4879 abu= 0.8179943
time= 160 ref. Vol= 189.4879 abu= 0.7704978
time= 161 ref. Vol= 21.03084 abu= 3.899036
time= 162 ref. Vol= 42.06168 abu= 3.61374
time= 163 ref. Vol= 42.06168 abu= 3.471093
time= 164 ref. Vol= 42.06168 abu= 3.233347
time= 165 ref. Vol= 13.48131 abu= 7.269324
time= 166 ref. Vol= 26.96262 abu= 6.750087
time= 167 ref. Vol= 26.96262 abu= 6.342114
time= 168 ref. Vol= 26.96262 abu= 6.156672
time= 169 ref. Vol= 10.29108 abu= 9.717158
time= 170 ref. Vol= 20.58215 abu= 9.134128
time= 171 ref. Vol= 20.58215 abu= 9.5714
time= 172 ref. Vol= 20.58215 abu= 9.085543
time= 173 ref. Vol= 7.73765 abu= 15.63782
time= 174 ref. Vol= 2.803497 abu= 36.02644
time= 175 ref. Vol= 5.606993 abu= 34.59965
time= 176 ref. Vol= 5.606993 abu= 32.99451
time= 177 ref. Vol= 5.606993 abu= 31.03268
time= 178 ref. Vol= 1.704253 abu= 62.7841
time= 179 ref. Vol= 3.408506 abu= 57.50319
time= 180 ref. Vol= 3.408506 abu= 54.86274
time= 181 ref. Vol= 3.408506 abu= 53.10244
time= 182 ref. Vol= 0.8876319 abu= 112.6593
time= 183 ref. Vol= 1.775264 abu= 108.1529
time= 184 ref. Vol= 0.6699109 abu= 155.2445
time= 185 ref. Vol= 1.339822 abu= 146.2881
time= 186 ref. Vol= 0.4481009 abu= 229.8589
time= 187 ref. Vol= 0.4481009 abu= 263.3335
time= 188 ref. Vol= 0.8962018 abu= 218.7007
time= 189 ref. Vol= 0.8962018 abu= 181.8787
time= 190 ref. Vol= 0.8962018 abu= 152.8674
time= 191 ref. Vol= 1.792404 abu= 101.5396
time= 192 ref. Vol= 1.792404 abu= 56.34892
time= 193 ref. Vol= 5.377211 abu= 30.12714
time= 194 ref. Vol= 5.377211 abu= 19.15491
time= 195 ref. Vol= 16.13163 abu= 10.78626
time= 196 ref. Vol= 16.13163 abu= 7.438801
time= 197 ref. Vol= 32.26327 abu= 4.21532
time= 198 ref. Vol= 64.52653 abu= 2.820545
time= 199 ref. Vol= 64.52653 abu= 1.54975
time= 200 ref. Vol= 193.5796 abu= 0.9453476
